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SOME CHARACTERISTICS OF A BOIUNG-WATER-COOLED 


VAPOR TRANSPORT FUEL PIN 
by James R. McConaghy, Jr. 

Lewis Research Center 

SUMMARY 

In this report the behavior of a boiling-water-cooled vapor transport fuel pin is in- 
vestigated. This type of pin is provided with a central void that permits fuel material to 
move from regions of high heat generation (and temperature) to regions of lower heat 
generation. This characteristic tends to remove hot spots and can be exploited to some 
advantage. A heat-transfer analysis is used to study pin characteristics at startup and 
during steady-state operation. A comparison with a conventional solid fuel pin of simi- 
lar design is made. 

The study investigates the fuel pin in the region of boiling burnout heat flux. Fuel 
pin and coolant channel dimensions are fixed. 

A comparison is made between vapor transport pins and conventional solid pins. 

For the particular case studied the average power density attainable in a vapor transport 
pin was 260 percent greater than in a conventional solid fuel pin. 

It is shown that under certain conditions a vapor transport fuel pin can be operated 
in the film boiling heat-transfer regime without the clad material reaching failure tem- 
peratures. This suggests that it may be possible to produce superheated steam from a 
properly designed single stage reactor core. 


INTRODUCTION 

The use of vapor transport fuel pins has been suggested (ref. 1) as a means of ex- 
tending the operating life of high power reactors. The vapor transport fuel pin is essen- 
tially a small pressure vessel with the fuel deposited on the inside surface. There is a 
relatively large central void. 

The pin is designed to operate such that the inside surface of the fuel is at a high 
enough temperature to permit vapor phase mass transfer of the fuel within the void. The 



movement of fuel within the void space tends to keep the inside fuel surface temperature 
constant and thus eliminate hot spots. In addition, the void space provides a storage 
volume for the gaseous fission products, thereby minimizing the rate of fission gas 
pressure buildup. A more detailed description of the vapor transport fuel -pin concept 
is available in reference 1. The major conclusions reached in reference 1 were that a 
vapor transport type of fuel pin would have longer life and produce higher fuel burnup 
than a conventional fuel pin under similar operating conditions. 

This report investigates the behavior of a single vapor transport fuel pin in a 
boiling-water -cooled reactor. Boiling-water reactors are of considerable practical 
interest because of their relative simplicity and because of the high pin surface heat 
fluxes attainable. In addition, boiling-water -cooled reactors are, at least to a degree, 
self-controlling (ref. 2). The design of an entire reactor system is not discussed here. 

Two computer codes are used to simulate operation of the fuel pin. The first code 
simulates pin startup with an initial uniform fuel distribution (this analysis is also used 
for the solid pin calculations) . The second code calculates pin characteristics at the 
final steady-state operating conditions after the fuel has redistributed to produce a uni- 
form internal fuel surface temperature. Figure 1 illustrates the initial and steady-state 
fuel distributions. Only radial heat conduction is considered within the pin. External 
heat transfer is to a flowing boiling-water coolant stream. 

This report presents a study of a 2.03-centimeter (0.800-in.) outside-diameter 
vapor transport fuel pin that is 122 centimeters (4 ft) long. The operation of this pin is 
studied through the region of the critical heat flux. Critical heat flux is defined herein 
as the point at which departure from nucleate boiling or "burnout" occurs. 

The vapor transport fuel pin is compared with a solid fuel pin under similar operat- 
ing conditions. The power distribution and power densities of the two types of pins are 
compared and discussed. 

This report does not fully investigate the behavior of the vapor transport fuel pin 
under all conditions. The main purpose here is to consider one particular case in order 
to illustrate some of the characteristics and advantages of designing a fuel pin to oper- 
ate as a vapor transport pin. Emphasis is placed on the heat-transfer aspects of the 
pin's operation. An analysis of the fuel pin lifetime as a function of deisgn parameters 
is available in a previous report by Rom, Finnegan, and McDonald (ref. 1). 


SYMBOLS 

A cross sectional area of coolant annulus 
c 

Cp heat capacity of coolant 


2 




Qw 

Q2 

^f-c 

^LB 


% 

R„ 


R„ 


R-, 


w 


heat capacity of coolant at fuel-pin surface temperature 
heat capacity of vapor phase coolant 

heat generated per unit thickness of fuel per thermal neutron 
thermal neutron attenuation coefficient 

temperature dependency coefficient of liquid phase coolant density 
constant of integration 

equivalent hydraulic diameter of coolant annulus 
fuel enrichment 

liquid phase convective heat-transfer coefficient 
vapor phase convective heat -transfer coefficient 

index used to identify axial increment in fuel-pin heat-transfer calculation 
thermal conductivity of fuel 
thermal conductivity of clad 
thermal conductivity of bulk coolant 

thermal conductivity of coolant at fuel-pin surface temperature 
fueled length of fuel pin 

total number of axial increments used in fuel-pin calculations 
power produced in axial increment Az of fuel pin 
coolant system pressure 

heat generation rate per unit volume of fuel at the exterior fuel surface 
heat generation rate per unit volume of fuel 
heat flux at the fuel-clad interface 
local boiling heat flux 

fully developed nucleate boiling heat flux at the value of T - T at which local 
boiling first occurs 

heat flux at the fuel-pin surface 

radius of coolant annulus 

radius of central void 

outside radius of fuel-pin clad 

inside radius of fuel-pin clad 
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r distance from fuel-pin centerline 

temperature of internal fuel surface 
T temperature of external clad surface 

A 

Tgat saturation temperature of coolant 

T, temperature of internal clad surface 

w 

T2 fuel temperature 

Tg clad temperature 

bulk coolant temperature 
V linear velocity of coolant 

z distance from coolant inlet end of fuel pin 

Az length of one axial increment used in heat -transfer calculation 
a volume fraction of vapor in coolant 
e effective emissivity 

X latent heat of vaporization of coolant 

Id viscosity of coolant 

Id viscosity of coolant at fuel -pin surface temperature 
p density of coolant 

Pj density of fuel 

density of liquid coolant at saturation 
p density of coolant at fuel-pin surface temperature 

O 

p^ density of vaporized coolant at saturation 
CT Stefan- Bo Itz man constant 

(p thermal neutron flux at fuel-clad interface 

unperturbed peak thermal neutron flux at fuel-clad interface 
X mass fraction of vapor in coolant 

Superscripts: 

(1) convective heat-transfer regime 

(2) transition region from convective heat transfer to fully developed nucleate boiling 

(3) fully developed nucleate boiling regime 

(4) film boiling regime 
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indicates an arithmetic average 


ANALYSIS 

The mathematical analysis of the fuel pin is developed in two parts. First, the fuel 
distribution is assumed to have uniform thickness for the entire pin length. Using this 
initial fuel distribution, a known axial neutron flux profile and known coolant inlet con- 
ditions, temperature and power profiles, are calculated. The first case predicts the 
maximum temperature and power level at system startup. 

The second case determines the required fuel load and fuel distribution necessary 
to produce a desired constant internal fuel surface temperature as a function of the peak 
heat flux. The program also has an option that will compute the internal temperature 
corresponding to a specified fuel load after fuel redistribution. The temperature and 
power distributions are then determined for this final operating condition of the pin which 
would exist after fuel redistribution. 

Both cases are steady-state solutions; one is the initial condition, and the other is 
the final condition. 


General Development of the Heat Conduction Equations 

The analysis is formulated with the following restrictions and assumptions: 

(1) End effects are negligible. 

(2) Fuel distribution and pin geometry are axisymmetric. 

(3) Heat conduction in the axial direction is negligible. 

(4) An average radial thermal conductivity for each material or region that is repre- 
sentative of the true temperature -dependent situation is used. This radial -average 
thermal conductivity is recalculated for each axial location. 

(5) The internal void space temperature at any given axial position is constant and 
equal to the internal fuel surface temperature, T^^, at that position. 

A schematic diagram of the pin model is given in figure 2. 

Based on restrictions (1) to (5) the heat transfer in the fuel shell is considered first. 
For a thermal conductivity that is constant in the radial direction at any given axial 
location and for conduction in the radial direction only, the heat conduction equation is. 



dr 


-rQgW 


( 1 ) 
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With the boundary condition 


r 


dTgCr) 

dr 


r=R. 


= 0 


( 2 ) 


integration gives 


dTgCr) 

dr 



rQgWdr 


Integrating again from r = R. to r R results is 

1 w 


w 





rQ2(r)drdr 


(3) 


(4) 


where T is the temperature at the fuel-clad interface. 

The heat conduction equation that describes the heat flow through the clad is 



(5) 


Integrating once results in 


r 


dTgCr) 

“dr 


= c 


( 6 ) 


where c is a constant. 

The heat flux at the fuel-clad interface is 


dTgCr) 

^f-c - dr 



dTgCr) 

“dr 



(7) 


Therefore, 


c 


-^w^f-c 




( 8 ) 
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and 


dTg(r) ^ -R^Qf.c 
dr rkg 


( 9 ) 


Integrating from r = to 


r = Rq yields 


T 


s 



( 10 ) 


The heat flux at the fuel-clad interface is simply the power generated in a 

length (increment) of fuel Az divided by the interfacial area or 


^f-c 


--■Z(J1_ 

27tR Az 
w 


if power is constant over Az. The power generation is given by 


( 11 ) 


P(J) = 27tAz 



rQ2(r)dr 


( 12 ) 


Equations (4) and (10) to (12) completely describe the radial heat transfer in the pin. 

A computer code used to calculate the characteristics of the pin at startup is de- 
scribed in the next section. 


Case 1: Description of Computer Code Used for Startup Calculations 

The fuel pin is broken up into N axial increments such that the values of the axial 
(z) coordinate index J run from J = 1 at z=0 to J=N at z = L (see fig. 2). The 
value of z associated with the increment is the midpoint of that increment. The 
J = 1 and J = N increments have a length of ^ Az. The value of N is chosen such that 
increasing the value has little effect (within ±1 percent) on the calculated pin operation. 

The heat generation Q2(i’) is a function of the thermal neutron flux <p{z). The flux 
profile is assumed to be a truncated sinusoid modified by the volume fraction of vapor 
a in the coolant channel such that 

(^) [o- 201 + 0- 070^ (r)] (13) 
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where is the peak flux that would reach the fuel if there were no voids in the coolant 
channel and where e is the enrichment. Equation (13) is assumed to be representative 
of a typical axial flux profile in a boiling water reactor. 

The heat generation rate at the fuel-clad interface Q^(z) is taken as 

Q^(z)=C20(z) (14) 

where C 2 is a function of the fuel type. For the purpose of this report the heat- 
generation rate is assumed to vary with radial position as 

-C„(R -r) 

where Cg is a thermal neutron flux attenuation coefficient. The use of a slab geometry 
equation simplifies integration and is valid for the thin fuel shell t 3 q>ical of a vapor 
transport fuel pin. 

The power generated by any fuel segment of length Az is then 


P = 27tAz 


Qw(2) 


/ 

*'R, 


w -C 3 (R -r) 
re ^ dr 


(16) 


This particular expression may be integrated analytically. The computer code is set up 
to evaluate the integral numerically so that any form of a power integral may be used. 

The pin is assumed to be surrounded by coolant in a channel of annular cross sec- 
tional area A The coolant properties are calculated on the basis of energy added to 
the coolant in the (J-1) increment. If 


^4(J-1) ^ '^sat 


(17) 


then the coolant is subcooled and 


-T , ^(J-1) 

4(J) - T4(J-1) p(j_i)Cp 


(18) 


P(J) - ^4^(J-1) 


(19) 


litliAJzll 

(J) P(j) 


(20) 


8 



If 


'^4(J-D 


> T 


sat 


( 21 ) 


and 


U 

^P(J) - ■^c’'^(l)^(l) [^p(^sat - '^4(1)) •*• 


( 22 ) 


J=1 


then the coolant is at the saturation temperature and is set equal to The 

coolant density and velocity are then given by 


(J) 


^(J-l)P(J-l) 




(23) 


V 


(J) 


^(J-l)P(J-l) 

^(J) 


(24) 


If 


u 

^P(J) >^cV(l)^(l) [^p(Tsat - '^4(1)) + 


J=1 


(25) 


then the coolant is superheated and 


4(J) 


r J 

■'^c^(l)P(l) f^p(^sat ■ '^4(1)) + 

^ J=1 


P(J) + ^cV(l)^(l)^P^ "^satf 


"^c^(l)^(l)^P. 


^( J )=^( J - 1)(-^1 


(26) 

(27) 
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must be absolute 


and is calculated as in equation (24) where and 

temperatures. 

The steam quality on a mass basis in the coolant channel at increment J is 



The void fraction of steam on a volume basis is 


^ , P{J) - 

Py Pl 


( 28 ) 


(29) 


Except for the first increment where T^j^ is read in, is first assumed to 

be equal to Then using equation (4) 


T = T. - 
w 




rQg (r)drdr 


(4) 


R R 


the temperature at the fuel-clad interface T can be calculated. Then 

AV 


^s 


w 



( 10 ) 


is used to calculate the fuel-pin surface temperature T . 

By using the calculated surface temperature and the calculated coolant temperature 
and physical properties (eqs. (17) to (29)), a heat-transfer rate to the coolant may be 
determined. 

This value is compared with the rate necessary to dissipate the energy generated by 

fission. If the two heat fluxes do not match within ±0. 1 percent, a new value of Ti(j) 

is assumed, and T /yv is recalculated. When the heat fluxes match within ±0. 1 percent, 

th 

the code continues to the (J+1) increment. 

Appropriate subroutines are used to calculate the heat-transfer rate to the coolant 
and to calculate temperature dependent physical properties. These are discussed in 
appendixes A and B. 
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The code is written to accept all data in SI units except pressure, which is express- 
ed in atmospheres absolute, and temperature, which is in °C. 


Case 2: Description of the Computer Code Used for 
the Fully Redistributed Fuel Calculations 

The indexing of the axial increments of the fuel pin is done in the same manner as 
for the startup calculation (case 1) . 

The thermal neutron flux distribution, heat generation rate, and power output are 
calculated as they were in the startup code. The appropriate relations are 


(cAs) o (l - f ) [»■ 201 »■ S’D" (r)] 

Q^(z) = C2<p(z) 

-CgCR -r) 

Q2(r,z) =: Q^(z)e 


P(j) = 27 tAz Q^(z) 


/ 


re dr 


(13) 

(14) 

(15) 

(16) 


Again, the pin is assumed to be surrounded by coolant in a channel of annular cross 
sectional area A . The coolant properties are calculated on the basis of energy added 

C 4.U 

to the coolant in the (J-1) increment. The method of calculation is shown in equations 
(17) to (29) in the previous section. 

This code is written to operate in two different modes. Mode one calculates the fuel 
loading necessary to produce a preselected void temperature. Mode two calculates the 
internal void temperature obtained from a preselected fuel loading. 

Mode one was used for most of the calculations in this report and is the basis for 
the following discussion. 

An internal temperature T^ is input. This temperature should be high enough to 
permit a reasonable rate of vapor transport in the void space (2600° C for UO2 fuel is 
used in this report) . 

Except for the first increment where Rj^j is read in, is first assumed to be 

equal to Then using equation (16), 
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P(j) - 27 tAz 


/ 

•'R. 


re or 


(16) 


i*h 

the power generated in the J increment is calculated. • The heat flux that must be 
transferred from the pin surface is then 


lo(J) 27 tAzR 


(30) 


The heat -transfer subroutine is then called (see appendix A) and the surface temperature 
necessary to obtain a heat flux of is calculated. 

The temperature of the inside surface of the fuel clad T^ can then be obtained by 




substituting ~ into equation (10), 





(31) 


By rearranging equation (4) the internal temperature T^ can then be obtained 




if'/ 




-c 




dr dr 


(32) 


i(J) ^i(J) 


If this internal temperature does not agree to within ±0. 1 percent with the desired tem- 
perature, a new value of Rj(j) is selected and T^^(j) is reculated. When the calculated 
and desired internal temperatures compare within ±0. 1 percent, the code continues to 
the (J+1)^^ increment. 

The fuel mass in each increment is calculated and summed to give the total mass 
of fuel necessary to produce the selected internal temperature. 

If the operation is in mode two, then the calculated fuel mass is compared with the 
desired fuel mass. If they do not compare within ±0. 1 percent a new value of the desired 
internal temperature is selected, and the entire mode one process is repeated. 

The temperature dependent physical property subroutine is discussed in appendix B. 

Again the code is written to accept all data in SI units except pressure, which is in 
atmospheres absolute, and temperature which is in °C. 
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DISCUSSION 


A study of a vapor transport fuel pin has been done using the computer codes de- 
scribed in the analysis section. The simplifying assumptions made in the development 
of the mathematical model should result in conservative calculations of pin behavior. 
Although refinements in the model could be made, the approach used here should not 
result in any gross errors. 

Curves are presented to illustrate the major characteristics of the vapor transport 
pin. A comparison between a vapor transport pin and a similar solid-fuel pin is made 
and the implications discussed. 


General Description of the System 

The study presented in this report is based on a 122 -centimeter (4-ft) long, 2.03- 
centimeter (0.8-in.) outside -diameter pin. The clad is 0.254-centimeter (0.1-in.) thick 

AISI 347 stainless steel. The fuel is uranium dioxide. 

0 

The coolant is water at a pressure of 6.9x10° newtons per square meter (68 atm) in 
an annular channel. The pin is in a vertical position with upward coolant flow. 

The unperturbed neutron flux at the fuel-clad interface is assumed to be a chopped 
sinusoid (axially) and symmetric in the angular direction. The actual flux at the 
fuel- clad interface is a function of the void fraction of steam in the coolant channel, 
as described by equation (13). In this equation, the flux is assumed to be proportional to 
the mass of water in the surrounding coolant flow channel. This model is intended to 
give an axial flux variation that is somewhat representative of the actual variation that 
would exist in a boiling water reactor. For any particular reactor, the actual axial flux 
variation should be used in the code. 

A reference case was selected as a basis for the calculations. The values of the 
major parameters for this case are as follows: 


Pin length, L, cm 

Outside clad radius, R^, cm 

Inside clad radius, R,„, cm 

w 

Radius of coolant annulus, R„, cm 

Coolant pressure, P^,N/m 

Coolant velocity, V, cm/sec 

Coolant inlet temperature, T^, °C 

Unperturbed peak neutron flux, c|)q( 0. 93/e), neutron/cm^sec 
Internal fuel surface temperature, Tj^, °C 


. . 122.0 
. . 1.015 
. . 0.762 
. . 1.602 
6.895X10® 
. . 305.0 
. . 275.0 
1.6X10^^ 
. 2600.0 
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The values of the constants used in the calculations are as follows; 


Heat generation rate per thermal neutron, Cg, cal/cm-m 1.0174x10"^^ 

Neutron flux attenuation coefficient, Cg, cm~^ 4.05 

Density variation coefficient of liquid phase coolant, C^, dimensionless 1.0 

Number of axial increments in computer calculation, N, dimensionless 100 

Density of fuel, pj, g/cm^ 10.96 


Effect of Neutron Flux 

The unperturbed peak neutron flux (p^{0.93/e) where e is the enrichment, is used 
as input data to the code as a means of changing the heat generation rate per unit volume 
of fuel. Figure 3 is a plot of the peak internal and external clad temperatures as a func- 
tion of 0^(0. 93/e). The calculation is for the reference case with fuel redistribution 
and a fuel loading sufficient to produce a constant internal fuel surface temperature of 

2600° C. The sharp discontinuity in the curves between <b (0.93/e) = 1.7x10^^ and 
13 ^ 

0^(0. 93/e) = 1. 8x10 neutrons per square centimeter per second indicates a transition 
of the coolant from fully developed nucleate boiling to film boiling ("boiling burnout 
point”). As neutron flux level increases, fuel loading must decrease to maintain con- 
stant fuel internal surface temperature (2600° C). However, the heat flux increases 
until film boiling burnout occurs. 


Initial Conditions 

Calculations were made to determine the initial temperature distribution for a fuel 

pin loaded with 1500 grams of uranium dioxide (the fuel load of the reference case), 

uniformly distributed on the pin walls as shown in figure 1. The pin was assumed to be 

1 3 

exposed to the same neutron flux level as the reference case (0 (0.93/e) = 1.6x10 

2 ° 
neutrons/(cm )(sec) unperturbed peak flux) . The resulting internal temperature profile 

is plotted in figure 4. The temperature spike at the center of the pin is due to film boil- 
ing of the coolant. Obviously the temperatures calculated would never be reached in 
reality because of redistribution of the fuel due to vapor transport. However, several 
important conclusions can be drawn from this initial-condition calculation. 

(1) If the fuel is not shaped initially to approximate the steady-state operating fuel 
profile, attention must be given to startup rates. 

(2) Internal temperatures high enough to cause vapor transport of the fuel (~2600° C) 
can be reached before boiling burnout conditions are reached at the pin outer surface. 
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(3) The method of calculation presented here should be extended to investigate the 
transient behavior of the fuel pin. 


Production of Superheated Steam 

The present analysis indicates that it is possible to design a system to produce 
superheated steam. One particular case is illustrated in figure 5. The data for this 
calculation are identical to the data for the reference case of the parametric study, with 
the exception of the coolant annulus radius which is reduced from 1. 602 to 1. 300 centi- 
meters. The reduced mass flow rate of coolant results in a higher void fraction of 
steam being generated. This flattens the neutron flux profile. Film boiling, due to the 
high steam void fraction, occurs at all flux levels in figure 5. 

By operating a reactor in this manner it may be possible to produce superheated 
steam in a single-stage system without exceeding clad temperature limitations. 


Comparison of the Vapor Transport Pin With a Similar Conventional Pin 

The power density produced by a conventional fuel pin almost completely filled with 
fuel (0. 02 mm central void) was compared with that of a similar pin operating as a vapor 
transport type of pin. Two constraints were put on the operation of the pins: (1) the 
centerline temperature must remain at or below 2600° C, and (2) the critical heat flux 
must not be exceeded. The neutron flux level is a dependent variable in this calculation. 
It assumes whatever values are required by one of the two preceding constraints. In 
practice, the first constraint limits the operation of the conventional pin, and the second 
constraint limits the operation of the vapor transport pin. 

The power density based on total pin volume attainable by both types of pins is plot- 
ted against outside pin diameter (constant clad thickness) in figure 6. The vapor trans- 
port pin permits much higher power densities than the conventional pin at pin diameters 
above 1.015 centimeters (0. 4 in.). In the example, both pins are filled with fuel below 
a diameter of 1. 015 centimeter (0. 4 in.), and the curves coincide. The power density 
drops off only slightly with increasing diameter for the vapor transport pin. This sug- 
gests that if vapor transport pins were used, fewer pins of larger diameter could be 
used to construct a reactor without a major penalty in power density loss. The resulting 
economic benefit could be significant. 

The higher power output of the vapor transport fuel pin is a result of two mechan- 
isms. One, the internal fuel surface temperature of the vapor transport fuel pin is kept 
at 2600° C for most of the pin length. Also, the fuel layer is much thinner than in the 
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conventional type pin. Therefore, because the insulating effect of the fuel is reduced, 
a higher surface heat flux without melting of the fuel is possible. Two, the heat flux 
profile for the vapor transport fuel pin is flatter in relation to the neutron flux profile 
than the heat flux profile for the conventional pin (fig. 7) . 

The vapor transport pin was shown in reference 1 to have other advantages. The 
central void of a vapor transport pin would extend the pin life significantly by slowing 
down the rate of pressure buildup due to the averaging out of the fission product gas 
accumulation. Also, the vapor transport of fuel would permit higher burnup fractions 
than possible in the conventional pin because of favorable replenishment of fuel from low 
reactivity worth zones to high reactivity worth zones. The dimensions and fuel loading 
of the example pins discussed in this report have not been optimized with respect of pin 
lifetime . 


CONCLUSIONS 

The performance of a vapor transport type of fuel pin in a boiling-water -cooled 
reactor has been studied using two computer codes. One code calculates pin startup 
(initial) conditions with uniform fuel distribution. The other code calculates final pin 
conditions at steady state with the fuel redistributed to produce an isothermal (constant) 
internal fuel surface temperature. 

The fuel pin analyzed in this report is 122 centimeters (4 ft) long with an outside 
diameter of 2.03 centimeters (0.8 in,). The uranium dioxide fuel is deposited on the 
interior walls of the pin cladding with a void space running axially for most of the length 
of the pin. The pin is operated such that the interior fuel surface is at a temperature of 
2600° C. Vapor phase mass transport within the central void maintains a constant tem- 
perature over the internal fuel surface independent of external variations in neutron flux 
or heat transfer to the coolant. 

The major conclusions of this study are 

1 . The vapor transport fuel pin can operate at a higher power density than a conven- 
tional pin. This effect is particularly noticeable at the larger pin diameters. Since a 
major percentage of reactor cost is in fuel pin fabrication this effect alone could result 
in a significant economic benefit. Other factors that affect power reactor costs are be- 
yond the scope of this report. 

2. Under certain conditions a vapor transport fuel pin can be operated above the 
critical heat flux without the clad material reaching thermal failure temperatures. This 
suggests that superheated steam could be produced from a properly designed single- 
stage reactor. 
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3. The startup rate for a vapor transport pin would be determined by the time re- 
quired for vapor transport to occur without burning out the pin. This startup time could 
be reduced by making the initial fuel distribution approximately the same as that for the 
final steady-state operating condition. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 16, 1970, 

122-28. 
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APPENDIX A 


BOILING HEAT TRANSFER CORRELATIONS USED IN THE FUEL PIN CALCULATIONS 


Case 1 : Heat Transfer Calculations for the Fuel Pin Startup Code 


The calculated fuel pin surface temperature (eq. (10)) is used to obtain a surface 
heat flux from a heat flux as a function of temperature driving force curve for boiling 
water. 

The curve is divided into the four regions which are described in the following sec- 
tions (ref. 3). Figure 8 illustrates the relation of the four regions to each other. 

Region 1. - This is the forced-convection regime where 


= hi (T^ - T^) 

The forced-convection heat -transfer coefficient hj is given by 


(Al) 


h 


1 



where the physical properties are those of the bulk coolant stream (ref. 4) . 
Region 2. - Region 2 is the partial nucleate boiling regime where 


(A2) 




1 +< 


^(3) 

.(i) 




(A3) 


This is an interpolation formula for the boiling curve in the transition region between 
forced convection heat transfer and fully developed nucleate boiling. 

The term is the fully developed nucleate boiling heat flux at the surface-coolant 
temperature difference at which local boiling first occurs. It can be obtained graphically 
as shown in figure 5-4 of reference 3. 

The onset of local boiling is given by the intersection of the forced convection curve 
(eq. (Al)) and the expression for local boiling heat flux q^^g where 


q^B = 1.1X10-4 




8 (Tg - T 


2.3/(14.7Pp)°-°^^^] 


sat'' 


(A4) 
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Note that equation (A 4) requires that pressure be in atmospheres and temperature be in 
°C to give a heat flux q^^g in kilowatts per square centimeter. 

Region 3. - This is the fully developed nucleate boiling regime where 


= (14. [l-S (Tg - Tgat)]^^ (^5) 

Equation (A5) requires that pressure be in atmospheres and that temperature be in °C 
to give in Btu /(hr) (ft^) kilowatts per square centimeter. 

Region 4 . - Region 4 is the film boiling regime where 


<!<•»' = (T^ - T,) + ec (T^ - 

Note that absolute temperature (K) must be used in equation (A6) . 
The heat-transfer coefficient h^ is calculated from 



(A6) 


(A7) 


All physical properties in equation (A7) are determined at the temperature of the pin 
surface. 

Between regions 3 and 4 there is another transition region between fully developed 
nucleate boiling and film boiling. Since information is not reliable for the region, it was 
omitted in the code. The transition to film boiling is treated as a discontinuity in the 
boiling curve which results in conservative calculation of surface temperature. 

The critical heat flux determines the point at which fully developed nucleate 

boiling changes to film boiling. For the purposes of this study a critical heat flux rela- 
tion was fabricated from available data to give values representative of experimental 
data in the range of interest. The relation is 

/V 

w = (A8) 

Equation (A8) is for use in this study only and should not be used as a design equation. 

In this study transition to film boiling also occurs if the steam fraction of coolant, 
on a volume basis, exceeds 0.80. 
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Case 2 : Fully Redistributed Fuel Code Heat-Transfer Calculation 


The heat flux as a function of temperature driving force curve is constructed as in 
case 1 using equations (Al) to (A8). In case 2, however, a surface temperature is ob- 
tained from a known value of the surface heat flux. 
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APPENDIX B 


TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY CORRELATIONS 


In the calculations presented in this report the thermal conductivity is allowed to 
vary with temperature in the axial direction. This appendix includes the correlation 
used to obtain thermal conductivities for the uranium dioxide fuel, the stainless -steel 
clad, and the coolant water as a function of these temperatures. 

UO2 Fuel. - Except for the first axial increment, the temperature used is the 
arithmetic mean of the inside clad temperature and the inside fuel temperature of the 
(J-1)^^ increment increment for J=l). The thermal conductivity correlations are 
(ref. 5) 

kg = -13.86x10"® T2 + 21,284x10"^ Tg ^ 667 K (Bl) 

kg = -3. 56x10"® Tg + 11. 442x10"^ 667 K < Tg < 1220 K (B2) 

kg = 6. 1702x10"® Tg > 1220 K (B3) 

The average fuel temperatures are all in K. The units on kg are (cal) (cm)/ 
(sec)(cm®)(°C). 

Stainless-steel 347 clad. - Except for the first axial increment, the temperature 

used is the arithmetic mean of the inside and outside clad temperature for the (J-1)*' 
th 

increment (J*"” increment for J=l). The thermal conductivity of the clad is (ref. 6) 


kg = 6. 75x10"® Tg + 33. 25x10 ® 


(B4) 


cy 

The units on k^ are cal/(sec)(cm j(°C/cm) for Tg given in K. 
Coolant water . - The temperature of the coolant is considered to be 
The thermal conductivity of the coolant water is (ref. 7) 


T4(j)in«C. 


k^ = 8.15x10"'^ T4 + 1. 515x10"® 
k^ = -1. 164x10"® T^ + 1. 805x10"® 
k^ = 4. 76x10"® T^ + 2. 702x10"® 


T^ < 150° C 

(B5) 

150° C < T^ < 250° C 

(B6) 

250° C < T4 - Tg^^ 

(B7) 
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1. 63X10-* 


The units on k. are (cal)(cm)/(sec)(cm^)(°C). 
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Coolant 

channel 


Figure 2. - Schematic diagram of fuel pin model used in 
mathematical analysis. 


Initial uniform Steady-state fuel 

fuel distribution distribution 

Figure 1. - Fuel distribution before irradiation 
and after the fuel has been fully redistributed. 
Data are that of the reference case. 
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Ratio of mass of fuel Maximum fuel pin surface temperature, 

in pin to total ca- 
pacity of pin 


Internal clad surface temperature 
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Figure 3. - Effect of neutron flux level. 



Axial position, z/L 

Figure 4. - Initial fuel pin temperature dis- 
tribution at startup for conditions of reference 
case. Fuel is uniformly distributed. 
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Figure 5. - Maximum surface temperature and total power out- 
put as functions of unperturbed thermal neutron flux. 

Fuel, uranium dioxide; coolant annulus radius, 1.300 centi- 
meters. (The coolant annulus radius is 0.302 centimeter 
less than that of the reference case. 


1 


Power generated, kW/m 



1.0 1.5 2.0 

Outside diameter of pin, cm 


Figure 6. - Comparison of vapor transport pin with conventional 
fuel element. Critical heat flux is not exceeded for either pin. 
Clad thickness, 1 inch. 
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Axial position, z/L 
lb) Conventional solid fuel pin. 

Figure 7. - Neutron flux at fuel surface and 
surface heat flux distributions for both 
vapor transport and conventional solid fuel 
pins. Pin diameter, 0. 8 inch (2.030 cm); 
reference case conditions. 
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Figure 8. - Boiling heat transfer curves for conditions at coolant Inlet end of refer- 
ence pin. 
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